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We present a class of observables which are suitable for determining the fidelity of a state to the 
multipartite Greenberger-Horne-Zeilinger (GHZ) state. Given an expectation value of an observable 
belonging to the class, we give a simple formula that gives a lower bound and an upper bound for 
the fidelity. Applying the formula to the GHZ-state preparation experiment by Pan et al. [Nature 
(London) 403, 515 (2000)], we show that the observed state lies outside of the class of biseparable 
mixed three-qubit states. We also show that for this class of operators, adopting the principle of 
minimum variance [Phys. Rev. A 60, 4338 (1999)] in the state estimation always results in the 
state with the minimum fidelity. 
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I. INTRODUCTION 

Recently, importance of entangled states of multipar- 
tite system has been realized not only as a fundamen- 
tal concept of quantum mechanics 0] but also as an es- 
sential resource for quantum information processing 0. 
Up to now, there have already been several experimen- 
tal reports on 3- and 4-particle entangled states 0, 0, @ 
One of the important measures to analyze how close the 
produced quantum state is to the desired maximally en- 
tangled state is the fidelity 0, i.e., the overlap with the 
desired entangled state. For bipartite systems and mul- 
tipartite systems, the fidelity is used as a criterion for 
nonseparability and distillability of the so-called Werner 
state or Werner-type state, i.e., a maximally entangled 
state mixed with the completely depolarized stateQ. In 
the classification of mixed three-qubit states it was 
shown that one of tripartite witnesses can be used to de- 
tect a state that does not belong to the biseparable class. 
This witness is given by 

W=il-P GHZ , (1) 

where Pghz is the projector onto a GHZ state. When 
Tr[W p] < 0, the state p lies outside of the biseparable 
class and has genuinely tripartite entanglement. Since 
Tr[W p] is written as 1/2 — / using the fidelity / to the 
GHZ state, the fidelity is useful to determine to which 
class a state belongs 9]. 

As can be seen from these examples, it is important 
to determine the fidelity / of experimentally produced 
states from the observed data. Systematic arguments 
about the possible fidelity values allowed by experimen- 
tal data for multipartite system will be helpful to the 
experimental realization of various applications in quan- 
tum information processing. The main purpose of this 
paper is to give a formula for the possible range of the 
fidelity value in the simplest case, i.e., the case where an 
expectation value of a single operator is given as exper- 
imental data. We present a class of observables which 



are suitable for determining the range of the fidelity of 
a state to the n-partite GHZ state |3> n ). The class is 
determined through the expansion of the projection op- 
erator | into the sum of direct products of Pauli 
operators for each party. Given an expectation value of 
an observable belonging to the class, the formula gives a 
lower bound and an upper bound for the fidelity. As an 
example, we analyze the GHZ-state preparation experi- 
ment by Pan et al. Q and show that the fidelity to the 
GHZ state is larger than 0.71. This indicates that the 
observed state does not belong to the biseparable class 
and has genuinely tripartite entanglement. 

In addition to the argument of what kind of states 
are possible under the constraints of experimental data, 
there also is a problem of determining which state is 
most likely under the constraints. Such problems of es- 
timation for bipartite system has been discussed along 
the maximum entropy principle |10|. Application of only 
this principle sometimes leads to an estimated state that 
possesses stronger entanglement than the minimum en- 
tanglement that is compatible with the measured data. 
Based on the additional assumption that the realization 
of a stronger entanglement is less realistic, Horodecki 
et al. introduced a new constraint, i.e., minimization 
of entanglement 11] in applying the maximum entropy 
principle. They thus obtained an estimated state that 
has the minimum entanglement. Rajagopal derived the 
same state with a different assumption together with the 
the maximum entropy principle, i.e., to minimize the 
variance of a Bell operator|l2j. Since then, much at- 
tention has been paid|l3| to this problem. Here we will 
show that, in multipartite systems, applying the min- 
imum variance principle to the operators belonging to 
the above class gives the states with the minimum fi- 
delity that is allowed by the constraints. This is a gener- 
alization of Rajagopal's results to multipartite systems, 
and reveals why and in what cases the minimum variance 
leads to small entanglement. 

This paper is organized as follows. In Sec. [HJ we take 
a GHZ state as the state of interest and present the class 
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of operators by decomposing the projector into the sum 
of operators forming a commutative group. In Sec. II I II 
we derive a simple formula that gives a lower bound and 
an upper bound for the fidelity. In Sec. IIVI we apply 
the formula to analyze the GHZ-state preparation exper- 
iment by Pan et al. Q, and show that the fidelity to 
the GHZ state is larger than 0.71. Sec.^deals with the 
state estimation problem based on the minimum variance 
principle. Sec. I VII concludes this paper. 

II. DECOMPOSITION OF PROJECTOR AND 
CLASS OF OBSERVABLES 

In order to discuss the fidelity, we have to specify one 
state of interest. We have referred to the state as the 
desired state, we take n-partite GHZ state as the 
desired state, which is defined as 

~ ^d+i; +2; • • • ; +») + l-i; -2; • • • ; -»»• (2) 

In the following discussion, we consider the fidelity of the 
state p to \$n), i-e., / := Tr[p |$ n )($ n |]. The projection 
operator |$ n )($„| can be expanded as 

1 / n n 

V=i 3=1 

n n \ 

3=1 3=1 / 

= (l/2 n )(O +O 1 + ---+O 2 n- 1 ), (3) 
where O p is defined by 

n 

3=1 

where the n-bit sequence bgbi 1 ' ' is the binary rep- 
resentation of p, b n — Y™—i bj, since the terms with odd 
parity for b\ ■ ■ ■ b n vanish in the above expansion. The 
superscript j of the Pauli operators denotes particle j. 
It is easy to see that O p O q = O p © 9 , where p © q is the 
bitwise XOR of p and q. Hence the set of 2 n operators 
{Op} forms a commutative group isomorphic to [Z-if 1 . 
We denote this commutative group as A„. The operator 
Oo is the identity operator for the 2™-dimensional space, 
and the other operators Oi, • • • , O^-i have two eigen- 
values, ±1. All elements of A n take |$ ra ) as an eigenstate 
with eigenvalue 1. 

For n = 2, the above expansion is explicitly written as 
follows, 

1*2) ($2! - (1/4) (I 1 / 2 + a\a\ + io\iol + a>*),(5) 

where I represents the identity operator for the 2- 
dimcnsional space. For n = 3, it is written as 

1*3) (*a| = (l/8)(% + Izz + O zIz + O zzI 

~~\~&xxx ~\~ @xyy ~t~ &yxy &yyx)i (6) 



where we have used the simplified notations as Om := 
7W, Izz := Wo*, O zIz := gPf O zzI := 

^z^z^ 1 Oxxx '■— < J x a x <J x'i ®xyy ■— O x %0 'yi&y , Oyxy - = 

ivl<rliVy, and O yyx := ia\ia\a 3 x . 

Now let us consider the problem of determining the 
fidelity by measuring the expectation value of an ob- 
servable. The most direct approach is, of course, to 
measure (|$„)($„|), which will be done by conducting 
2™ — 1 different correlation measurements to determine 
(Oj) (j — 1, 2, • • • , 2 n — 1). Our interest here is how we 
can deduce the information about the fidelity from ob- 
servables that can be measured by much smaller number 
of correlation measurements. In what follows, we give a 
formula to derive an inequality for the fidelity from an 
expectation value of an observable belonging to the class 
C n defined as follows. 

Class C n : An observable A belongs to C n if 
and only if A is a linear combination of operators 
Qii 0,2, • • • , Qm € A„ with positive coefficients, and the 
set {Qi, Q,2, ■ ■ • , Qm} forms a system of generators for 
A„. 

The class is determined through A„, and hence deter- 
mined by the desired state |$ n )- 

The minimum cardinal number of systems of genera- 
tors for A„ is n. To see this, suppose there are only l(< n) 
generators. In this case, however, they can generate at 
most J2i=oiCi = 2 l {< 2") kinds of elements of A„ due 
to the property of the Pauli matrices. An example of a 
system of generators for A„ is 

n-1 

{O X,X, - ■ : X i O z ,I,I,- ,z, 0/, z , ,z, ■ ■ ■ , OlJ t ... ,z,z}> (7) 

where X}X> ... lX := alal ■ ■ -a x , 0/, z ,/,..., z := 
I 1 all 3 ■ ■ ■ <7™, and so on[lJ|- This system of gener- 
ators for A„, indeed, generates all 2" elements of A„ 
with the help of Eq. JI). 

As for A2, examples of a system of generators are 

r 1 2 1 21 rl2 12 • 1 * 21 Ct 12 -1-21 
Wx<Tx, <*x a z}l Wx<, Vz^z, IVylVy}, {* , °x°xi ™y™y} 

and so on. As for A3, examples of a system of genera- 
tors are }(D X yy, Oyxyi Oyyx}i \®xyyi Oyxyi @Izz}: \&xxxi 
Oxyyi @yxy> Oyyx} and SO On. 

III. INEQUALITY FOR FIDELITY 

In this section, we derive an inequality under the con- 
dition that the expectation value of an operator A € C n 
is specified. First, we will show that if the given expecta- 
tion value is the maximum value, the state must be \& n )- 
For that, we use the following lemma. 

Lemma. Let X, Y be operators taking eigenvalues ±1, 
and [X,Y] = 0. Then, 

1-\(X)-(Y)\>(XY)>(X) + (Y)-1. (8) 

This is directly proven by inequalities: 

((1- X)(l-Y)) > and ((1 ± X)(l =p Y)) > 0. (9) 
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In the following, we consider A — aiQi + a 2 Q 2 + • • ■ + 
ctmQm (2™ > m > n, ati > 0, V i). We assume that 
A G C n . We assume that the maximum expectation value 
(A) = M is given, where M := X)t=i a i ^ s the maximum 
eigenvalue of A. Because all coefficients at are positive, 
(A) = M implies (Qi) = (Q 2 ) =■■■ = (Q m ) = 1. Since 
Qj — I for any j and the set {Qi, Qi-, - ■■ , Qm} forms a 
system of generators for A„, any element Q of A„ can be 
written as 



(10) 



with 1 < 0i < 02 < • ■ ■ < 01 < m. We show that (Q) = 1 
for all I as follows. When I = 1, then (Q) = 1 holds. 
Suppose (Q) = 1 holds when I = k, i.e., (n»=i QPi) = 1- 
Then with the help of the lemma, (rL=i Qft) = 1 holds. 
Let X be [[ i=1 Qf3i and Y be Qp k+ll respectively. The 
lemma leads that 1 < QPi) — 1- Hence (Qi) = 

(Q2) = • • • = (Q m ) = 1 means the expectation values of 
all elements of A„ are one. This means that the fidelity 
is 1 with the help of Eq. Therefore we obtain the 
following: 

Proposition 1: Let M be the largest eigenvalue of A S 
C n . If Tr[p A] = M, then p = |$„)($„|. 

This implies that the largest eigenvalue M is not de- 
generate. This is crucial point in deriving an inequal- 
ity for the fidelity. Let us write the eigenvalues of A as 
M, r 2 , r 3 , • • • , and r s (s < 2"), where M > r 2 > r 3 > 
• ■ • > r s > —M. In this notation, when some eigenval- 
ues of A are degenerate, then s < 2" holds, and when 
all the eigenvalues of A are not degenerate, then s = 2™ 
holds. Since proposition 1 implies that |$ n ) is the only 
eigenstate for the largest eigenvalue M, we can generally 
expand (A) := Tr[p .A] as 



(11) 



i=2 



where / = Tr[p |$ n )($„|] is the fidelity to |$„), and 
5i(> 0) satisfy / + 2*_ 2 ?i = 1- Using this relation to 
eliminate 92, we have 



(A) = M/ + (l-/-^ ?i )r 2 + X>r i 

?' — 3 i—3 

s 

- (M-r 2 )/ + r 2 -^g i (r 2 -r i ), (12) 



i=3 



We thus obtain 



M-r 2 



(13) 



The equality of the relation (|13() holds when $^f_ 3 ?j = 0. 
We can also derive an inequality that gives an upper 



bound of the fidelity by eliminating q s , namely, 

s-l s-l 

(A) = Mf + J2<Hri + (l-f-J2^ r ° 

i=2 i=2 
s-l 

= (M-r s )f + J2<li( r i-r s )+r s , (14) 



and 



{A) r, 
M-r s 



(15) 



The equality of the relation l|15|) holds when J2i=2 1i = ^ - 
We therefore obtain the following proposition. 

Proposition 2: Let M, r 2 , and r s be the largest, the 
second-largest, and the smallest eigenvalue of A € C n , 
respectively. When (A) := Tr[p A] is given, the fidelity 
/ := Tr[p I |] is bounded as 



(A) - r 2 ^ (A) - r s 
M-r 2 ~ M- r s 



(16) 



IV. APPLICATION TO EXPERIMENTAL DATA 

We analyze the experimental data by Pan et al.Q. In 
this experiment they obtained four expectation values of 
three-photon polarization correlations, 



(xyy) 



(yxy) 

(xxx) ~ 0.74. 



(yyx) ~ 0.70, 



(17) 



These experimental data is obtained by the post selec- 
tion, i.e., picking up only the events with each of the 
three detectors registering a photocount. If the detectors 
used in the experiments were ideal ones, the post-selected 
state would be contained in a 2 3 -dimensional subspace, 
which can be identified with a tripartite system of three 
qubits. Then the above observed values could be con- 
sidered to give the expectation values (O xyy ), (O yxy ), 
{Oyyx), and {Oxxx)- In the real experiment, however, 
the detectors are not ideal, namely, they cannot distin- 
guish a single photon from more than one photons and 
they have a limited quantum efficiency and dark count- 
ing. Due to these imperfection together with the non- 
ideal photon source, the post-selected state also contains 
contributions outside of the 2 3 -dimensional subspace, in 
which two photons or no photons enter the same detector. 
However, the superfluous contributions can be neglected 
as compared to the statistical uncertainty (a few %) of 
the observed values as follows 0. The contribution of no- 
photon events are due to dark counting, but the rate of 
the dark counts is low enough [3| to be able to neglect 
the effect. The contribution of more than one photon en- 
tering a detector passes the post selection only if another 
detector have a dark count, or more than two photon 
pairs are created in the parametric downconversion. The 
former case is negligible due to the low dark count rate, 
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and the latter is also negligible since the probability per 
pulse to create n-photon pairs is of the order of about 
1CP 4 ". We can thus assume that the post-selected state 
is related to polarization of three photons and approxi- 
mately supports 2 3 -dimensional Hilbcrt space. Hence we 
obtain 

(O xyy ) ~ (O yxy ) ~ (O yyx ) ~ 0.70, 

(O xxx ) ~ 0.74. (18) 

The limited quality of the polarization optics just be- 
fore the detectors may make the visibility lower, which 
will make the estimated fidelity smaller. Hence we use 
these experimental expectation values for restricting the 
fidelity from below. 

Clearly, each observable of (|18|l does not belong to C3 . 
Therefore we take the summation of these expectation 
values. We then obtain 

{A) ~ 2.84, (19) 

where A = O xyy + O yxy + O yyx + O xxx . Apparently, A <E 
C3 holds, and eigenvalues of A are 4, 0, and —4, where 
is a degenerate eigenvalue. With the help of Eq. (|TB|l . 
where the parameters are set to M = 4 and r% = 0, 
we can state that the observed state in this experiment 
has the fidelity to a GHZ state larger than or equal to 
0.71 [l^. Because the value is lager than 1/2, with the 
help of Eq. 0}, these experimental data, indeed, ensure 
the observed state does not belong to the biseparable 
class and has genuinely tripartite entanglement. 

V. RELATION BETWEEN VARIANCE AND 
FIDELITY 

In this section, we show that if we apply the minimum 
variance principle for estimating the states from an ex- 
perimentally obtained expectation value of an operator 
belonging to C n , the estimated fidelity becomes the min- 
imum value that is allowed by the constraints. This is a 
generalization of Rajagopal's results 01 to multipartite 
systems. 

Let us consider Rajagopal's case, i.e., bipartite sys- 
tem in which the variance of Bell-CHSH operator B = 
\/~2(a x a x + <jI<Jz) is made minimal. 

Suppose that an expectation value (A) is given, where 

Case 1 A = v x <J x , 

Case 2 A= a\a 2 x + a\a\. (20) 

Note that, for Case 1, the operator A does not belong 
to C 2 , whereas, for Case 2, the operator A belongs to C 2 , 
and A = B/V2. 

For Case 2, combining the results by Refs.^l| and 
[l2|. it is shown that if variance ((AA) 2 ) is made mini- 
mal, the calculated fidelity to |<£> 2 ) takes minimal when 
the expectation value (A) (> 0) is given. In this way, 
the minimum entangled state was derived from Jaynes 



principle. For Case 1, the given (A) determines a range 
< / < ^ A \ +1 for the possible value of /, and the mini- 
mization condition for ((A^4) 2 ) puts no further condition 
on this range. The allowed fidelity value is thus unsettled 
and can be any value in the region, except for the case 
that (A) = -1 (-> / = 0). (Remember -1 < (A) < 1 
for Case 1). This also means that, for Case 1, we cannot 
derive minimum fidelity (i.e., zero) from the minimum 
variance principle. 

Next we consider several examples for tripartite sys- 
tem. There are eight GHZ states, which are written as 

|Vi(8)) := (l/V2)(|+i; +2; +3) ± |-i; -2; -3)), 
1-02(7)) := (l/V2)(|-i; + 2 ; +3) ± |+i; -2; - 3 », 
1-03(6)) := (l/\/2)(|+i; - 2 ; +3) ± hi! +2; -3)), 
104(5)) := (l/V2)(l+i; +2; -3) ± |-i; -2; +3».(21) 

The eight probabilities of observation of these GHZ states 
for the state p are defined as 

Pi := P |V<>, (t = l,2,-»,8). (22) 

We can see that is equal to The fidelity / is 

then identical to pi . 

Suppose that an expectation value (A) is given, where 
we consider 

Case 1 A (D xyy ~\~ & yxy . 

Case 2 A = O xyy + O yxy + Ozzii 

Case 3 A — (D xyy 4~ (D yxy -t- (D yyx . (23) 

Note that, for Case 1 and 2, the operator A does not 
belong to C3, whereas, for Case 3, the operator A belongs 
to C s . 

For Case 3, when the variance ((A^4) 2 ) is made min- 
imal, we see, later, that we can derive the minimum fi- 
delity. For Case 1 and 2, however, the minimum vari- 
ance principle does not work completely. In Case 1 and 
2, the given (A) determines range < / < M±H and 
< / < 4 1 , respectively for the possible value of 
/. For Case 2, the minimization condition for ((A.4) 2 ) 
puts no further condition on this range |l6j. For Case 
1, the minimization of ((A„4) 2 ) makes / max smaller, i.e., 
the allowed fidelity value is for —2 < (A) < and 
< / < ^# for < (A) < 2. (Remember -2 < (A) < 2 
for Case 1). This means that, for Case 2, we cannot de- 
rive the minimum fidelity (i.e., zero) from the minimum 
variance principle and for Case 1, we cannot derive the 
minimum fidelity for < (A) < 2. 

Now we calculate the fidelity for Case 3 from the min- 
imum variance principle. If we write operator A in the 
matrix form using the GHZ basis, the diagonal elements 
become 3, —1, — 1, —1, 1, 1, 1, and —3, and no off- 
diagonal element appears. This means that the measured 
value for A can take four values 3, —1, 1, and —3, where 
— 1 and 1 are degenerate eigenvalues. Using notations 
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Pa ■= P2 + P3 + Pi, P/3 ■■= P5 + P6 + P7, and the rela- 
tion f = Pi, the probability that the measured value for 
A takes 3, 1, — 1, or —3 is expressed as {f,P/3,Pa,Ps}- 
We can calculate this for the three cases corresponding 
to 3 > {A) > 1, 1 > (A) > -1, and -1 > (A) > -3, 
as follows. When (A) lies between 3 and 1, the mea- 
sured value for A can take only 3 or 1 but not —1 or —3 
to attain the minimum variance of its distribution. The 
minimization of ((A„4) 2 ) thus leads to the distribution 

{/) Pp, Pa, Pa} to be j 3 ~^ ,0,0|. Similarly, 

for 1 > (A) > — 1 case, the distribution is calculated to 
±±M izM }, and for 





■(A) 



-1 > {A} > -3 case to 



be {o, 2 , 2 
be{0,0,2±M,^V 

by the minimum variance principle, is then summarized 
as 



j. The derived fidelity in Case 3, 



/ = 



(A) -I 
2 





3 > (-4) > 1 
1 > (A) > -3. 



(24) 



It is easy to show that this is equal to the minimum of 
the fidelity values with the help of Eq. i|16|) . where the 
parameters are set to be M = 3 and r 2 = 1. It is thus 
concluded that, for Case 3, the derived fidelity to |$3) by 
the minimum variance principle leads to the minimum of 
the possible fidelity values allowed by the expectation 
value. 

We generalize the argument for n-partite system. We 
consider A — aiQi + U2Q2 + • ■ • + a m Q m G C n . Suppose 
that an expectation value (.4) is given. In the following 
we calculate the fidelity / from the minimum variance 
principle in this case. We write the eigenvalues of A as 
M> r i, ?"3, • • • , and r s , where M > r 2 > r 3 > ■ ■ ■ > r s > 
—M. The probabilities for observing these eigenvalues, 
M, r 2 , 7-3, • • • , and r s are denoted as /, 52, 93, • • • , and 
q s , respectively. Similarly to the discussion as to Case 
3 for tripartite system, if (A) lies between M and r2, 
the minimization of ((A.A) 2 ) leads to Yll=3 1i = 0) which 



(A)-r 2 M-{A) 



means {/, q 2 , <&,♦•• ,q s } = j M _ r2 , M _ r2 
The derived fidelity is then summarized as 



,0. 



(A) r 2 

f= ) M-r 2 



M > (A) > r 2 
T2 > (A) > r s . 



,0}. 



(25) 



Eq. I|24() is a special case of Eq. (|25|l where the parameters 
are set to be M = 3, r 2 = 1 and r s — —3. We can see 
that Eq. 125|l gives the minimum fidelity with the help of 
Eq. 17] . Hence we have the following result: 

Proposition 3: When (A) := Tr[p .4] is given, where 
A € C n , the derived fidelity to from the minimum 
variance principle is the minimum of the possible fidelity 
values allowed by the expectation value. 

VI. CONCLUSION 



In conclusion, we have analyzed the possible fidelity 
values that are compatible with an expectation value of 
a single operator as experimental data. We have defined 
the desired maximally entangled state and formulated 
one class that is related to a decomposition of the pro- 
jector onto the desired state. We have made use of the 
commutative group theory to formulate a class of observ- 
ables. When an expectation value of an operator that 
belongs to the class is given, we can derive an inequal- 
ity that gives a lower bound and an upper bound of the 
fidelity values that are compatible with the expectation 
value. With the help of the inequality, we have analyzed 
the experimental data by Pan et a/.|j]. The data ensure 
the observed state does not belong to the biseparable 
class and has genuinely tripartite entanglement. Finally, 
we have also analyzed the calculated fidelity from the 
minimum variance principle. 
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